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1. Introduction 

The CKM matrix element \V c b\ is a key quantity in Standard-Model tests. It normalizes the legs 
of the unitarity triangle. The dominant uncertainty in | V c t | comes from theoretical determinations 

of the hadronic form factors for B — > civ + The exclusive processes B — > D£v and B —> D*£v 

[1, 2, 3] can be studied in lattice gauge theory. Here we report on results for B — > Dtv. Lattice 
calculations at zero recoil typically have the smallest errors. However, because of the phase space 
suppression near zero recoil in B — > D£v, experimental errors are largest there. Thus, we aim to 
work at nonzero recoil where the combined experimental and theoretical error is minimized. This 
work updates our previous report [4] with all ensembles now included in the analysis. 

The differential decay rate dF(B — > D£v) /dq 1 is, for mj <C min(Mj,g 2 ), proportional to |/ + | 2 
for i = e, /J,, where for q = pg — Pd, 



(D(p D )\y»\B(p B ))=f + (q 2 



Ml-M 2 



q 2 



+ fo(q 2 ) B 2 q ■ (1.1) 



Here = by^c is the b — ^ c vector current and / + and fo are the vector and scalar form factors, 
respectively. The alternative form factors h + and /j_ are convenient: 

(D(P °jlgg ( ''- )> =/. + W(v + v-)^/,-(w)(v-,/r. (..2, 

where v = pb/Mb and v' = po/Mp. They are related to / + and fo through 

U(q 2 ) = ~[(i+r)M^-(i-^-ML 



fo(q 



2\ 



w + 1 w — 1 

h + (yv)-- /j_(w) 



(1.3) 



J+r " v ' 1 
where r = M D /M B and q 2 = M\ +Mp - 2wM B M D or w = v ■ V . 

2. Lattice-QCD calculation 

We are carrying out calculations on 14 lattice ensembles generated in the presence of 2 + 1 
flavors of improved (asqtad) staggered sea quarks [7] with light-quark masses and lattice spacings 
shown in Fig. 1. 

In the B meson rest frame for any recoil D-momentum p we can obtain h + and /?_ from matrix 
elements of the current starting from ratios of lattice matrix elements R + and R , and x/, where 

R + (p) ^ <D(p)|V 4 |B(0)> 

*_(p) " < D(P)|V1|B(0) > 



<D(p)MS(0)> 



(P(p)|W|P(0)) 

^ (p) ^ WTO) i2A) 

w(p) = [l+x f (p) 2 ]/[l-x f (p) 2 ] 
h + (w) =/?+(p)[l-x / (p)/?_(p)] 

Mw) =*+(p)[i-Mp)/*/(p)]. 
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Figure 1: Left panel: parameters of the 14 lattice ensembles in this study. Plotted are values of the light to 
strange sea quark mass ratio m^/m/, vs. the approximate lattice spacing in fm. Right panel: example from 
the a = 0.06 fm, m^/m/, = 0.15 ensemble with T = 24,25. Upper curve: the zero-recoil double ratio, middle 
curve: the ratio R + (p,t)/R + (0,t) for the IS smeared D-meson interpolator, and bottom curve: the local 
interpolator. Red points are included in the fit (p — 0.15). 



At zero recoil, we can also use the double ratio of Hashimoto et al [5] : 

(D(0)\V l \B(0))(B(0)\V l \D(0)) 
1 +WI (D(0)|V 4 |D(0)> (B{Q)\V 4 \B(Q)) ' K ' 

The continuum and lattice currents are matched through "V^ = Z v ^V^ b . We use a 
mostly nonperturbative method [5], writing 



Z v? b =Pv^ Z v*Zv* b , (2-3) 
and determine p v ii from one-loop lattice perturbation theory. 

cb 

In addition to the two-point functions, we need matrix elements (F(p)|V' I |X(0)) for X,Y G 
{B,D}. They are constructed from naive light spectator quark propagators and clover heavy quark 
propagators in the Fermilab interpretation [6] as shown in Fig. 2. Valence bottom and charm quark 
masses were tuned to the "kinetic" B s and D s masses. The mass of the naive light spectator quark 
is set equal to that of the light sea quark. 

For interpolating operators Ox, we compute two-point and three-point functions 



C 2 ^*(p,0 = ( O x (0)O x (t)) , (2.4) 
C 3 ^(p;0;,r) = (ol{0)V^{t)O x (T)) . (2.5) 



We use both point and IS smeared interpolating operators for the D meson and 15 smeared inter- 
polating operators for the B. 
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Figure 2: Valence quark line diagram for B — > D. 



2.1 Two-point and three-point correlator fits 

We obtain the lattice form factors via a two-step procedure. First, we fit the B- and D-meson 
two-point correlators to obtain the energies and overlap factors. Then we use these determinations 
as constraints (with Bayesian priors) in the three -point fits. For illustration, we show the reduction 
of the three-point and two-point functions to obtain R+(p) = (D(p)|V 4 |B(0)). We include excited 
B and D contributions indicated with a prime but not both together: 

4f^°(p,0 = v ^^<D(p)|V 4 |B(0)>- 7 = r v^(0) 

+ ^^)f^(D\p)\V 4 \B(0)) e -^^^{0) (2.6) 
\/2E D i \/2m B 

. e - E Dt -m B ,(T-t) 

+ ^ZMp)^==(D(p)\V 4 \B'(p)) v _ V^(0), 

or 

Cjf B ^ D (p,f) =C (p)(D(p)\V 4 \B(0))e- EDt e-"^ T -^ [l + C 1 {p) e - AEDt + C 2 (p)e ( - t - T ^ AmB ] , 

(2.7) 

where Cb(p), AEp = Ejy — Ed, and Am B = m B i — m B come from fits to two-point correlators. Terms 
oscillating as (— )' (not shown) are introduced by the naive light quark. We suppress their contri- 
butions by averaging over T, T + 1 and t,t + l, as introduced in [1]. 

Putting information from three- and two-point functions together, we get 



clT B ^ D (vM ED -' nB)t+(mB - mD)Tl2 j z D (0)E D 
'cl ! !- D ^ D (0,t)cl P 4 '^ B (0,t) i Z D{v)m D 



Ra 



l+s l (p)e- AEDt +s 2 (p)e^ A ' nB 



(2.8) 



The zero-recoil form factor h+(0) = ^+(0) can be calculated very accurately from the double ratio. 
A good strategy is to use it to normalize the nonzero recoil values: 

1T7^ = ^M^P(Smt)+A(p)exp(-AE D t)+B(p)exp(Am B t), (2.9) 

where 8m = 0, AEp = E D i —Ed, and Am B = m B > — m B are constrained by fits to two-point functions. 

We do a simultaneous fit to three three-point functions, as illustrated in the right panel of Fig. 1 
and determine R + , R , and Xf for each momentum (recoil parameter w), from which we determine 
h + and 
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Figure 3: Form factors h + (left) and /i_ (right) as a function of the recoil parameter w. The curves in each 
panel are a result of a simultaneous fit to the chiral-continuum expressions of Eq. (2.10) with p = 0.68 for 
the h + fit and p = 0.41 for the /i_ fit. 



2.2 Chiral-continuum extrapolation and q 2 parameterization 

The resulting form factors h + and /j_ are shown in Fig. 3. We fit them to the expressions 



h + (a,me,w) 



h-(a,m£,w) 



1 -p + (w- 1) +k + (w- 1) H ^ hci i+ m^ + c a)+ fl! + c a>w& a (w-1) 

(2.10) 



mf. 



g 2 * 

— ^02 lo g S l-loop( A Z » w > 



X 



in. 



pl{w- \)+k-(w- l) 2 + ci -m? + Ca.-a 2 + c a , w ,-a 2 (w - I) 



for light spectator quark mass m^, lattice spacing a, and w = v • V . For the one-loop chiral logs we 
use a staggered fermion version of Chow and Wise [9]. Thus, these fit functions contain the correct 
next-to-leading-order chiral perturbation theory expressions, including staggered discretization ef- 
fects [10]. As can be seen, the dependence of h + on a and m^/ni/, is quite mild. We expect h- 
to have larger discretization effects than h + because of different HQET power counting. This is 
consistent with what we see in the data. For | V c b\, the contribution coming from /j_ over the entire 
kinematic range is small, so the larger errors in /j_ don't increase the overall error much. These 
features with 14 ensembles are consistent with our previous findings with four ensembles [10]. 

To compare the lattice and experimental form factors we need to extrapolate to larger w (equiv- 
alently q 2 ). We do this using the z-expansion of Boyd, Grinstein and Lebed [11], which provides a 
model-independent parameterization of the q 2 dependence of / + and /b- This expansion builds in 
constraints from analyticity and unitarity. It is based on the conformal map 



z(w) 



Vl+w-v^ 
VI+W + V2 ' 



(2.11) 



5 



Semileptonic B to D decays at nonzero recoil 



Carleton DeTar 



B->D zexpansion x 2 /dof=5.9/3 
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Figure 4: Left: form factors /+ and /o parameterized by the z expansion (p = 0.12). Right: comparison 
with experimental results from the Babar collaboration [12]. The red line gives our result. The red dotted 
lines show only statistical errors. The boxed region appears to have the smallest combined error. 



which maps the physical region w € [1, 1.59] to z £ [0,0.0644]. It pushes poles and branch cuts far 
away at \z\ ~ 1. Form factors are then parameterized as 

^ CO 

m =™£°"' z "- <2 - 12) 

where P/(z) are the Blaschke factors and 0/ are the "outer functions". The latter are chosen to 
simplify the unitarity bound: 

£k„| 2 <l- (2-13) 

n 

In practice, we need only the first few coefficients in the expansion. We also impose the kinematic 
constraint /+ = fo at q 2 = or z ~ 0.0644. 

To implement the z expansion, we start from the value of / + and fo at the physical point, as 
determined from the chiral/continuum fit. We choose four w values, w = 1.00, 1.04, 1.10, and 1.16, 
and use the corresponding form factor values to determine the coefficients a,-,o, and a/2. These, 
then, are used to parameterize the form factors over the full kinematic range, as shown in the left 
panel of Fig. 4. 

We compare our result with experimental measurements from the Babar collaboration [12] in 
Fig. 4. For present purposes we take \V c b\ from B — > D*£v at zero recoil [3]. 



3. Future plans 

To complete the analysis, we need to apply small corrections resulting from adjusting the 
charm and bottom quark masses to their tuned values, implement the full current renomalization, 
and compile a complete error budget. 
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